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ABSTRACT 



Very high energy gamma-ray flares from the Crab nebular detected by AGILE and 
Fermi satellites challenge our understanding of the pulsar wind nebulae. The short 
duration of the flares, only few days, is particularly puzzling since it is much shorter 
than the dynamical times scale of the nebular. 

In this work we investigate analytically and via numerical simulations the electro- 
magnetic signatures expected from the large amplitude low frequency magnetosonic 
waves generated within the Crab nebular which induce the corrugation perturbations 
of the termination shock. As a result, the oblique termination shock produces time- 
dependent, mildly relativistic post-shock flow. Using the relativistic MHD version of 
the RIEMANN code, we simulate the interaction of the termination shock with down- 
stream perturbations. We demonstrate that mild Doppler boosting of the synchrotron 
emission in the post-shock flow can produce bright, short time scale flares. 



The Crab gamma-ray flare: constraints from observations 



The detection of gamma-ray flares from the Crab nebular by AGILE and Fermi satellites is 



one of the most astounding recent discoveries in high energy astrophysics Tavani et al. (2011); 



Abdo et al. (2011). The flares typically continue for several days/weeks. No changes in the pulsed 



emission of the Crab pulsar or a glitch were detected. No variability at has been reported by other 
satellites: INTEGRAL reported no detection of the flare the 20-400keV window (ATel#: 2856) 
and Swift/BAT did not see any significant variability during the gamma-ray flare in the 14-150keV 
range (ATel#: 2893). Swift also reported no evidence for active AGN near the Crab, suggesting 
that the Crab itself is responsible for the flare (ATel^^: 2868). The prevailing conclusion from 
the observations of flares is that the flares are associated with the nebular (and not the neutron 



- 2 - 



star) and are mostly likely due to the highest energy synchrotron emitting electrons. Thus, the 
flares reflect the instantaneous injection/emission properties of the nebular and are not expected 
to produce a noticeable change in the IC component above ~ 1 GeV (Bednarek & Idee 2010). 



The two most surprising properties of the flares are intermittency and short time-scale vari- 
ability. For large intervals of time the gamma-ray emission form the Crab is nearly constant, with 
large swings in emissivity taking place only sporadically. This suggests a rare, hard to predict 
event that has disproportionately high-impact, a "Black Swan" event, in the nebular as being the 
cause of the flare. The flare duration, only a few days, is two orders of magnitude smaller than the 
dynamical time-scale of the nebular and several times smaller than the light crossing time of the 
termination shock. This presents major problem in the interpretation of the flares as variations in 
the structure of the termination shock. 

In present paper we consider the effects of the downstream turbulence on the properties of 
the termination shock. The nebular behaves as a resonant cavity capable of sustaining oscillations 
over a range of frequencies and wavelengths. Non-linear interactions of these waves could give rise 
to turbulence. A three or four a fluctuation in the turbulence could give rise to a strong wave 
that interacts with the termination shock, giving rise to an intermittent flare. As we demonstrate 
analytically and via numerical simulations, for particular combination of the viewing angle, overall 
shock obliquity and the amplitude of the perturbing wave, the intensity of synchrotron emission in 
the shocked plasma can indeed experience very short spikes, on time scales much smaller than the 
period. 



Komissarov &: Lyutikov (2011) suggested that flares come from the so called inner knot, a 



Doppler-boosted emission from the high- velocity flow downstream of the oblique termination shocl<[^ 
of the pulsar wind ( Komissarov &: Lyubarsky||2004 ) . The higher resolution simulations of the Crab 
Nebula discovered strong variability of the termination shock, involving dramatic changes of the 
shock shape and inclination (Camus et al. 2009). Thus, both observations of the morphological 



features of the Crab nebular (Hester 2008), as well as numerical simulations (Bucciantini & Del 



ZannapOOG Camus et al.||2009 ) show dynamical behavior of times scales of months to years. This 
discovery suggests that the gamma-ray variability may be related to the changes in the character- 
istics of the Doppler beaming associated with this structural variability of the termination shock. 
In this paper we study the interaction of a fast magneto-sonic wave with a relativistic MHD shock. 



Another possible scenario has been suggested by Uzdensky et al. (2011 ) who posit a very rapid 



reconnection event. We, however, point out that for reconnection to be rapid, it must take place 



within a turbulent environment (Kowal et al. 2009). Without picking a speciflc mechanism, we 



point out that any intermittent event with a bursty energetic release will trigger strong magneto- 
sonic waves that interact with the termination shock. Owing to the indeterminacy, that interaction 



^ We use the terms normal and oblique to indicate the relative direction of the shock normal and the fluid velocity, 
and the term perpendicular for the relative direction of the shock normal and the magnetic field. 
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will be brief. 



2. Model set-up 



In this paper we consider the interaction of strong magnetosonic waves generated within the 
bulk of the nebular with the termination shock. In Appendix |A] we consider a general problem of 
relativistic oblique perpendicular magnetosonic shocks. The results of these calculations are used in 



^2.1 to give simple estimates of the post-shock flow velocity and the conditions required to produce 
narrow high brightness peaks. 

In ^we describe relativistic MHD simulations of the interaction of strong magnetosonic waves 
generated within the bulk of the nebular with the termination shock. Unlike the large scale simu- 
lations of the nebular by ( Komissarov &: Lyubarsky||2004 , Camus et al.||2009 ) we zoom- in onto the 
small scale details of the variability of the post-shock flow. As the underlying steady-state flow we 
chose the highly oblique part of the termination shock, corresponding to the inner knot, Fig. [l]( see 



also Fig. 1 of Komissarov Sz Lyutikov (2011)). The choice of the highly oblique part of the initial 



shock is important. In these parts the unperturbed post-shock flow is already mildly relativistic, 
with the post-shock bulk Lorentz factor of the order of ~ l/(7r/2 — 9). Then even mild variations 
of the post-shock parameters are expected to produce large variations of the Doppler factor. 



2.1. Theoretical outline 



The PWN can be considered a resonator with internal frequencies of the order of the fast 
mode travel time trough the nebular. The resonator will support normal modes with the wave 
length ~ 1/L, the size of the nebular and frequencies ~ c/L ~ months-years. As normal modes 
propagating in the shocked flow fall onto the termination shock, they will be reflected back into the 
post-shock medium (the reflection angle is not equal the incoming angle, the reflected amplitude 
may be damped or amplified). During the reflection, the waves distort the shape of the termination 
shock, which becomes corrugated, see Fig. [2] 

Studying the corrugation perturbations of shock is an important topic in modern fluid dynamics 



(Landau & Lifshitz 1959 Fowles & Houwing 1984). Barring the spontaneous emission of waves 



by the shocl^ the shock corrugation can be calculated as a reflection of the sound waves by the 
shock back into downstream medium. What we are interested in here is the non-resonant response 
of a shock to a sound wave impinging from the shocked plasma side. That is, we are looking not 
for normal modes of shock oscillations, but for a response of the shock to arbitrary, non- resonant 
perturbation. 



^ Under certain conditions a shock can spontaneously emit a sound wave 
requires non-trivial equation of state: in polytropic gas this does not occur 



(D'yakov 


1954 


Kontorovich 


1957 


Landau & Lifshitz||1959 1. 
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Fig. 1. — Interpretation of the inner knot as a highly obhque part of the termination shock (Komis- 



sarov Sz Lyutikov 2011). 



Since the reflected waves intensity depends on the relative phase of the incoming waves and 
the corrugation waves, the corresponding perturbations of the downstream medium will include an 
entropy wave and one cannot use the barotropic equation of state; this complicates the problem 



considerably (Courant &; Priedrichs 1948, §73). Still, the salient features of the interaction can 



be derived in the limit of small amplitude of corrugation, in which case the entropy wave is weak 



(Courant Sz Friedrichs 1948, §73). In the isentropic limit, the shock acts as a partially reflecting 



surface (Landau Sz Lifshitz 1959 , §91), with the angle of reflection not equal the incidence angle. The 
incident and the reflected wave then form an interference pattern in the shocked fluid. Propagation 
of the shock through this interference pattern then induces weaker shock corrugations. 

As a simpler estimate we next consider the post-shock flow from a given corrugated shape of 
the shock, without calculating the dependence of the corrugation on the amplitude of the perturbing 
wave and neglecting the the entropy wave. General conditions at relativistic oblique perpendicular 




Fig. 2. — Geometry of the model. Initially the termination shock is an oblique planar shock, 
so that the incoming velocity makes angle 9 with the shock normal. Magnetosonic waves coming 
from the post-shock flow induce corrugations of the shock surface (solid line). Locally, the shock 
is oblique, making and angle (f) with respect to the upstream velocity. The post-shock flow makes 
an angle x with respect to the initial velocity (the post-shock flow is not orthogonal to the shock 
plane) . 

MHD shocks are considered in Appendix [A| 

Let us next assume that for the unperturbed oblique shock the angle of attack is 9, see Fig. [2] 
In addition, a wave propagating from the downstream induces shock corrugations, so that in the 
frame of the shock surface is located at y = 6y cos kx, where 6y is the amplitude of the corrugation 
and k is the wave number of corrugation waves propagating in x direction; the upstream flow is 
along y direction. The angle of attack at point x is then 

cos 9 + A sin kx sin 9 

cos m = , ( 1 

v/l + (AsinA;x)2 

where A = k6y is dimensionless amplitude of corrugations. Using the oblique shock conditions 
( |A9l|AT3l ) we can then calculate the post-shock velocity and the deflection angle, see Fig. [s] The 



maximum post-shock velocity is reached at wave phase tt/2 and equals 



/ 9(A2 + 1) 
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This value can be much larger than the initial unperturbed post-shock momentum, A = 0, ti2,o = 
^9/(8 sin^ 9) — 1. For mild amplitudes of corrugation, A < 1, the post-shock velocity is especially 
high for highly oblique shocks, 9 ^ A. 

These estimates demonstrate that even for fluid oblique shocks, which have substantially non- 
relativistic post-shock velocity for normal shocks, oblique corrugated shocks can have at specific 
phases very high post-shock velocities. Flows with higher magnetization have even higher post- 
shock velocities. Fig. |4j 



0.4 
0.2 
>< 0.0 
-0.2 
-0.4 







\ A=l /I 




A=3 





20 



15 



10 



1 2 3 4 5 6 
phase 








3 4 
phase 



Fig. 3. — Deflection angle x (Isft panel) and post-shock momentum ui (right panel) as a function of 
the phase of the perturbing wave for fluid shocks cr = for different wave amplitudes A = 0.1, 1, 3. 
The initial obliquity angle is 9 = 75°. For sufficiently high amplitudes a given defiection angle can 
be achieved up to four times per period. 



Simulations 



In ^3.1 we describe the numerical set up and the parameter range we investigated in our 



numerical simulations. In ^ 3.2 we show the results of our simulations and discuss their implications. 



3.1. Numerical Methods and Simulation Parameters 



The simulations were run on a 480x240 zone mesh using the relativistic version of the RIE- 
MANN code for astrophysical fluid dynamics ( Balsara|l998 ). While the code contains higher order 



effects (Balsara et al. 2009), in this application we used a second order ADER scheme along with 



a HLL Riemann solver and a WENO reconstruction. The upper and lower y boundaries were 
periodic, the right (pre-shock) x boundary were inflow, and the left (post-shock) boundary were 
set to produce the inflowing magneto-sonic waves. 
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Fig. 4. — Post-shock momentum ui as a function of the phase for different magnetization parameter 
cr = 0, 0.1, 1 and wave amphtudes A = 1. 

In the relativistic version of the RIEMANN code we evolve the conserved variables with: 



dtU + diP = 



where 





p 




p 


u = 




p = 


w;. 




u 







and 



S = ph-f'^v + E X B 
U = ph^^ - P + ^{E^ + B^) 



(3) 
(4) 

(5) 
(6) 
(7) 
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We use a F-law equation of state which gives 

r p 

We evolve the induction equation as 

dtB' + e''^'''djEk = (9) 

with ideal MHD approximation E = —v x B and the magnetic field co-located at the zone bound- 
aries in order to preserve the diB^ = divergence free constraint. 

Here n is the Eulerian density (n=7yo), p is the Lagrangian density, v is the three- velocity, S is 
the three-momentum, 7 is the lorentz factor, E and B are the electric and magnetic field vectors, 
P is the pressure, h is the enthalpy, F is the polytropic index (F = | in all our runs), and U is the 
energy density. All vector quantities are defined in the Eulerian frame, i.e. the rest frame of the 
shock. 6j is the usual Kronecker delta and e*'-'''^ is the three dimensional Levi-Civita symbol. 

This application studies the long-term interaction of waves with a stationary shock. We, 
therefore, require shocks to stay perfectly stationary on a computational mesh for long intervals 
of time when they are unperturbed. The theory of stationary relativistic MHD shocks is only 
reasonably useful because analytically exact solutions for such shocks are only obtained in the limit 
of vanishing pre-shock pressure. Using the exact solutions calculated using the jump conditions 
in Appendix |A] as an approximate starting point, we evolved the solution with finite but small 
pre-shock pressures on a very long one dimensional computational mesh. Finite pressure and 
discretization errors do make this shock move slowly with respect to the mesh. By running these 
one dimensional simulations for a very long time, we were able to quantify the speed with which the 
shock drifts on the long mesh. We then lorenz boosted the fluid variables into the rest frame of the 
shock. The resultant densities, pressures, velocities and magnetic fields in the pre- and post-shock 
regions for the stationary shocks have been tabulated in Table [T| 

The resulting numerically stable post-shock values were used to calculate the eigenvector for 



the fast right-going magneto-sonic wave using the analytical procedure from(Balsara 2001). The 
resulting wave was scaled so the density perturbation was a fraction of the unperturbed post-shock 
density, and then allowed to propagate into the shock at a 45° angle to the shock normal. We 
explored a range of angles of attack {(p) to test the effect of the obliqueness of the shock, and 
magnetization (a) to examine the effect of the magnetic field strength in the pulsar wind. We also 
varied the amplitude ( Sp/pi) of the perturbing magneto-sonic waves. 

Table [T] describes the primitive variables we used to set up the shocks. In all cases the frame has 
been rotated so that the shock is in the y-z plane, with the magnetic field in the positive z direction. 
The odd numbered runs had perturbation strengths of Sp/pi = 0.3 while the even numbered runs 
used Sp/pi = 0.7. Table [2] gives the parameters defined in Appendix [A| corresponding to these 
choices of primitive variables and compares the defiection angles and compression ratios calculated 
using the numerical code and the analytical expressions. 
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Table 1: Run parameters used to set up the un-perturbed shocks for our simulations, p is the 
density in the Lagrangian frame of the fluid, Yx is the fluid velocity normal to the shock, Vy is the 
fluid velocity along the shock, and is the magnetic field perpendicular to the shock normal and 
the velocity vector. The velocities and magnetic field are all given in the rest frame of the shock. 



Run 


Po 




My 


Po 




Pi 




Pi 


Bzi 


1,2 


0.1118 


-0.8452 


0.5286 


0.001 


0.2878 


3.8050 


-0.2542 


9.100 


0.951 


3,4 


0.1047 


-0.7007 


0.7057 


0.001 


0.9363 


1.720 


-0.2674 


1.948 


2.466 


5,6 


0.1190 


-0.3482 


0.9270 


0.001 


1.1254 


0.6680 


-0.1538 


0.152 


2.5610 


7,8 


0.1144 


-0.9831 


0.1698 


0.001 


0.2946 


5.011 


-0.3054 


16.1850 


0.9510 


9,10 


0.1059 


-0.9523 


0.283 


0.001 


0.9468 


1.9960 


-0.3860 


2.8140 


2.3380 


11,12 


0.1488 


-0.8937 


0.4150 


0.001 


1.0683 


1.1446 


-0.4373 


0.724 


2.1810 



Table 2: Variables describing the geometry of the un-perturbed shocks. <j) is the angle of attack, x 
is the deflection angle, a is the magnetization of the pre-shock plasma, r/ is the compression ratio 
between the post- and pre-shock fluids, and 70,1 are the Lorentz factors for the pre- and post-shock 
plasmas respectively. The table also compares the deflection angles Xnum and compression ratios 
^nwmcalculated using the numerical code and the analytical expressions, Xan and rjan-, calculated 



Run 





~X.num 


Xan 


a 


Vnum 


Van 


70 


71 


1,2 


57.9776" 


32.2950° 


32.1166° 


0.0046 


32.2458 


33.850 


12.6806 


1.2347 


3,4 


44.7963" 


24.0439° 


24.0094° 


0.0922 


16.4279 


16.359 


9.5316 


1.5241 


5,6 


20.5872" 


11.1670° 


11.052" 


0.2068 


5.6134 


5.474 


7.1745 


2.9233 


7,8 


80.2006" 


19.2743° 


19.324" 


0.0036 


43.8024 


44.1785 


14.6140 


1.0673 


9,10 


73.4494" 


19.6968° 


19.7816° 


0.1103 


18.8480 


19.042 


8.7586 


1.1390 


11,12 


65.0917" 


18.5929° 


19.9364° 


0.2230 


7.6922 


10.1557 


5.8646 


1.2534 



3.2. Simulation Results and Analysis 

Our simulations all show initially strong corrugations when the fast magneto-sonic wave first 
impinges on the shock. Figure [5] shows the density for our Run 5 shortly after the incoming wave 
has hit the shock, clearly showing the corrugation of the shock front. Larger perturbation am- 
plitude produces more significant and persistent corrugation of the shocks. In extreme cases the 
perturbations could even drive the shock front off of the computational domain. These pertur- 
bations effectively constitute an extra ram-pressure term, which increases the effective post-shock 
pressure, moving the shock to the right. In most situations, however, the perturbations are gentle 
and the shocks either show a slow, secular rightward motion or none at all. 

As time progresses the reflected waves coming off the shock interfere with the incoming waves, 
and the corrugation of the shock front begins to damp away. This damping is due to the homoge- 
nization of the effective ram pressure term as many wave fronts hit the shock. Once the simulation 
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Fig. 5. — Density for run 5 at early time, shortly after the wave hits the shock. The corrugation 
wave in the shock front is especially pronounced at these early times. 
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Fig. 6. — y-velocity (transverse velocity) for Run 12 at later time. The reflected waves are inter- 
acting with the incoming waves and creating this fish-scale pattern. 



has stabilized we see a prominent fish-scale pattern formed by the incoming and reflected waves, 
with the largest variations in the fluid variables occurring at the intersections of the waves. Figure 
[6] provides an example of this, showing the transverse velocity in Run 12, after the reflected waves 
have had time to fully cross the post-shock fluid. 

In all our runs we see the strongest fluctuations in the flow variables, and correspondingly the 
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strongest variation in the observed intensity, at the early times when the shock is most deformed 
by the impinging magneto-sonic waves. After the shock front has had time to regain equihbrium 
the variation in the Lorentz factor is smaller, and we do not see the brief spikes in intensity from 
Doppler boosting. 



4. Modeling the Gamma-ray Emission 

In this section we discuss the factors responsible for the observed high energy gamma-rays in 
the instantaneous approximation. Later, in ^we discuss the effect of time-of- flight considerations 
on the observed intensity. 

The highest energy leptons emitting synchrotron radiation at ~ 100 MeV have a very short 
cooling time scale, of the order of tcooi ~ 10^ 7^ sec, where 7 is the bulk factor of the post-shock 
flow and the estimate is done for the inner wisp located at ~ 10^^ cm. Thus, we expect that only 
a fairly narrow region behind the shock contributes to the very high energy emission. In addition, 
the observed radiation is highly modulated by the relativistic beaming effects. One expects that 
intrinsic variations of proper emissivity have only a marginal influence on the observed flux: it is 
mostly dominated by Doppler boosting. 

As a simple prescription that takes into account the radiative decay of shock particles and the 
Doppler boosting we integrate along the line of sight the quantity 5^ (here 5 = 1/(7(1-/3 cos 9i)) is 
the Doppler factor, 9i is the instantaneous angle between line of sight and the fluids velocity) for a 
fixed downstream distance from the shock. This is a simplified prescription that captures the two 
salient features of the post-shock flow: fast radiative decays and the dominant effects of relativistic 
Doppler boosting. 

Figure [7] shows a rough estimate for the intrinsic synchrotron emissivity in Run 12, calculated 
using the approximation IcxP*B", with a=1.5 as in Jun & Jones (1999); Balsara et al. (2001). As 



with the other variables, we see the highest intensity in the regions around the intersection of the 
incoming and reflected waves. While synchrotron radiation is expected to produce the high energy 
gamma-rays, the effect of the line of sight dependent Doppler boosting should be the dominant 
term determining the observed intensity. 

Figure [8] shows the Lorentz factor for Run 6, which represents our most oblique parameters 
{(p = 20.59°) with strongest perturbations {5p/pi= 0.7). One can see that the perturbations in 
the post-shock flow (with unperturbed Lorentz factor 71= 2.9 ) are large, ranging from 2.09 to 
3.96. This large variation in the post-shock Lorentz factor can lead to strongly enhanced observed 
intensity at some preferred lines of sight angles due to Doppler boosting, see Fig. [9j 

In the relative intensity profiles calculated using this simple radiation model described above, 
we see generally broad profiles for most line of sight angles, with a few angles producing sharp 
peaks from small sections of the corrugation wave (Fig. [9]). These peaks could be as high as twice 
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Fig. 7. — Synchrotron emissivity (IocP*B°) for Run 12 at early time with a=1.5, calculated in the 
fluids rest frame. The effect of the variation in this proper emissivity is expected to be negligible 
compared to the effects of the angle dependent Doppler boosting. 
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Fig. 8. — Lorentz factor for Run 6 at early time. The strong variation (a factor of 2) in the post- 
shock Lorentz factor leads to strong variation in the Doppler boosting along some lines of sight, 
which leads to spikes in the intensity. 

the baseline intensity. Since the peaks are related to small segments of the shock front, they have a 
small characteristic time scale. These flares are most pronounced at early times of the wave-shock 
interaction, when the shock is strongly corrugated. 
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Fig. 9. — Calculated intensity as a function of the phase of the corrugation wave for Run 12 at 
early time, with the line of sight angle at -45° to the shock normal. This profile assumes a constant 
rest frame emissivity, and only reflects the effect of the LoS dependent Doppler boosting ( IccS'^ ) 

5. Time-of-flight effects: Cherenkov-type interference 

The main drawback of the simplified radiation modeling described above is that it does not 
take into account the time-of-flight effects. One expects that neglecting the time-of-flight effects 
smooths out the intensity variations. We expect that the full temporal analysis will produce even 
sharper peaks: there is a special combination of the incoming wave direction and the amplitude of 
the corrugation such that the emission from the particular point of the corrugation wave add up 
coherently. Let us demonstrate this for a simple case of initially normal fluid shock. 

Consider a wave falling from downstream onto ? relativistically strong normal shock at an angle 
a to the shock normal. In the moving fluid sound waves become dispersive. For an unperturbed 
fluid moving relativistically with velocity vi = 1/3 in the medium with sound speed Cg = l/VS, 
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the corresponding dispersion relation (Anile 1989, Eq. 10.18) becomes 

\/50 — 2 cos 2a — 3 cos a 



k 



13 



The phase speed of corrugation perturbations along the shock normal is then 

VSO — 2 cos 2a — 3 cos a 



Vph 



13 sin a 



(10) 



(11) 



The corrugations become luminal for incidence angles smaller than a = 0.32. For smaller incidence 
angles the corrugation wave is superluminal. If a post-shock flow at some phase satisfies the 
Cherenkov condition 

sinx = c/vph (12) 



the waves emitted by this phase will add up constructively. Equation (12) together with the 
expression for the phase speed (11) and the post shock flow angle (A13) with (p = Tr/2 deflnes the 
condition on the incidence angle of a sound wave coming from the shocked fluid, such that the 
resulting corrugations of the shock surface produce the post-shock flow which remains in phase 
with the corrugation wave. In this case the corrugation pattern acts as a real particle, emitting 
light at the Cherenkov condition. 



6. Discussion 



In this work we considered emission from a relativistic corrugated shocks. We performed 
relativistic simulations of the corrugations of the perpendicular MHD shocks induced by strong 
magnetosonic waves impinging from the downstream medium. Our main conclusions is that long 
wavelength perturbations with the relative amplitude of ~ 0.5 — 0.7 of the downstream density 
can result in sharp intensity variations, with duration of the peak much smaller than the wave 
period. Sharp intensity variations can have amplitudes of an order of magnitude in flux, with large 
modulation index, of the order of 0.5. 

Variable, mildly relativistic post-shock flow may also explain (or, in fact, be required in stochas- 
tic shock acceleration models) the very high break energy observed in the Crab flares. As argued 



by Lyutikov (2010) (see also |de Jager et al. 1996), there is an acceleration model-independent 
upper limit on the frequency of synchrotron emission by radiation reaction-limited acceleration of 
electrons: 



h 



mc 



100 MeV. 



(13) 



The break energy observed by Fermi satellite in Crab's quiescent state, ~ 100 MeV (Fermi Collab 



oration &: Pulsar Timing Consortium||2009 ), approaches this limit, while the break energies during 
flares, ~ 500 MeV (Tavani et al. 2011 Abdo et al. 2011), exceed it. If emission is generated in 



a relativistically moving plasma, the maximum observed synchrotron frequency is boosted by the 
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Doppler factor. Note also, that previously EGRET data did indicate a moderate level of the cut-off 



energy variability (de Jager et al. 1996). 



Another important feature of our model is that in the proposed model the time scale of the 
flare is set not by the radiative decay time of the emitting particles, but by the overall dynamics 
of the termination shock. Thus, the flare time scale cannot be used to estimate the magnetic field 
within the emission region. 

We would like to thank Sergey Komissarov for suggesting the possibility discussed in this 
paper. 
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A. Oblique relativistic MHD shocks 



We derive jump conditions at oblique relativistic perpendicular MHD shocks, when the mag- 
netic field is orthogonal to flow velocity which is not aligned with the shock normal. We investigate 
the post-shock parameters as functions of the upstream plasma magnetization a, Lorentz factor 71 
and the angle of attack cj). Previously, Webb et al. (1987) investigated oblique relativistic MHD 



shocks in the quasi-parallel case, when the planes of fluid velocity and of the magnetic field coin- 



cide. Konigl ( 1980, ) considered relativistic fluid oblique shocks. Here we consider relativistic oblique 
perpendicular magnetohydrodynamic shocks, so that the magnetic fleld is in the plane of the shock 



and is orthogonal to the fluid velocity. This appendix corrects an error in Lyutikov (2004). 



Let the stream lines make an angle 



an angle x with the initial velocity. Fig. 10 Normal shock corresponds to 



with the shock normal and let the post-shock flow make 

= tt/2. 



Keeping only x-dependence, the equations of perpendicular relativistic MHD read {e.g., Anile 
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y A Shock 
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Fig. 10. — Geometry of the flow 



1989) 



dx (nux) = 

dx [p + + {w + B^)ul) = 
5, {{w + S2)7n,) = 
dx (uB) = 



(Al) 



where n is proper density, w = e + p = n + T/{T — l)p is fluid enthalpy, e is internal energy, B is 
rest-frame magnetic fields, n is a four- velocity. Note that in this Appendix we use the rest-frame 



magnetic field and density as primitive variables, while the numerical simulations in ^3.1 use the 
quantities in the lab frame. 



Oblique shock conditions can be obtained from normal shocks and an additional condition that 
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Konigl 


1980 


Anile 


1989 



noMo sin = niui sin((j) — x) 
Bquq sin (j) = Biui sin((/) — x) 

po + Bll2 + (w;o + Bl)ul sin^ = pi + B\I2 + {w^ + BI)u\ - x) 

(■Wo + ^0)701*0 sin = (it;i + B\)-ixUi sin((/) - x) 

/3o cos = /3i cos(0 — x) (A2) 

where subscripts refer to upstream (0) and downstream (1) regions and /3 is three- velocity. Below 
we assume the relativistic equations of state with F = 4/3. 

The problem is to express the post-shock parameters {j)\^n\^B\,u\ and x) in terms of the 
upstream parameters {po,nQ, Bq,uo and (p). As we demonstrate below, the corresponding equation 
take especially simple form is expressed in terms of the compression ratio rj = ni/no, which in turn 
depends on the upstream parameters. By ideal condition Bi = ijBq. Compression ratio can be 
related to the turning angle x 

rj' = {l + ul sin^ - 70^ sin^c^ - x)) (A3) 

sm (</> - x) 



In term of rj, 



. {{pi-Po) + {Bf-Bi)/2)rj^ 



. 2,. , {{p,-po) + {Bl-Bl)/2) r,^ 
7]^{Wq +B^) - (wi +Bf) 
2 Uq sin^ (f) if cos^ 'Z' + To sin^ 

Ml 



if sin^((/) — x) ^(1 + ^0 ^^^^ ' 

sin^ 



s{T?{(t) - x) = T[ . 2 . (A4) 



r/^ cos^ cf) + sin^ ^ 



Introducing magnetization parameter a = Bq/wq, the jump conditions (A2) can be resolved 
for post-shock pressure pi and enthalpy wi as functions of rj: 



Wocr{r]'^-l) 2 ■ 2 ,/i , ^^1 /l + WoSin\ 
Pi = PO 7, ^ ^0^0 sm 0(1 + 0-) 1 - W „ " ■ 2 



= W + (1 + ct)W 2 ^ °2 ■ 2^1 (A5) 

In case of cold plasma, po = the compression ratio r] can then be simply related to the inflow 
parameters: 



iHl-^nJ ? : + \ + . + 4„g(, + 1) s,n' ^ ■ + 1 1 I - „ - 2. ^ (AO) 
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Let us next simplify the above relations in the high compression limit, ry — t- oo. (The strong 



compression limit becomes inapplicable for sin cp < l/(no\/r+~fT).) Eq. (A6) gives 



16(1 + cr) sin^ 



[I -9a) sin^ 



8(72 + 2 + 161 + 0- + 5) (7 + VI60-2 + 16(7 + 1 + 1 



1 sin 



(7 <1 
(7 > 1 

(A7) 



The turning angle is given by 

(1 — Z) cos sin ( 



tanx 



cos2 (h + Z sin^ . 



(72 + 2 Vl6(72 + 16(7 + 1 + 5 cr + Vl6(72 + l6cr + 1 + 1 



1/3 < 2: < 1 (A8) 



^ 8(72 + 2 (^Vl6(72 + 16(7 + 1 + 13j (7 + Vl6(72 + 16(7 + 1 + 17 

Maximum deflection angle is reached at cos (p = sj Z jiX + Z and equals tan Xmax = (1 — ^)/ (22) 



l/2(tan(/) - cot0), Fig. 11 
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Fig. 11. — Maximum deflection angle Xmax as a function of the magnetization parameters a in the 

limit Uq OO. Xraax = 7r/6 for (7 = 0, Xmax ~ 1/(4(7) for (7 > 1. 
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The post-shock four-velocity 



2 • 2 

Ui sm : 



(9-8 sin^ 6) + 



8cj2 + 2 ( Vl6o-2 + l6cj-hl-F- 4 ) 0- 716^2 + iGcj + 1 - 1 



16(a + l) 



(A9) 



12 



(for normal shock, </) = tt/2, Eq. ( |A9| ) reproduces Eq. (4.11) of Kennel & Coroniti 1984), see Fig. 
In the high a limit, uj sin^ (j) ^ a + (9/8 - sin^ (f)) -\- 0{l/a). 
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Fig. 12. — Dependence of uisini;^ (the post-shock four-velocity times angle of attack) the angle of 
attack (p for different magnetization parameters o" = 0, 1, 10 (bottom to top). 



Using the post-shock fast velocity (Harris 1957) 



4 



Bf + 4/3pi 



4pi + m 



(AlO) 
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the post-shock Mach number is 



_vl _ 3{Z + l){{Z^-l) sin^0 + l) 



7Z - 1 

1 I 8 cos^ (t> , ( 2 56 cos^ (f>\ ^ ^ ^ i 

3 + + 3 — 3 r^' ^ ^ 1 



1 + 



The flow becomes supersonic for 



(see Fig. 13). 



cos^ 0-2/3 



0- > 1 







(All) 



(A12) 




Fig. 13. — Dependence of the deflection angle x on the angle of attack (j) for different magnetization 
parameters a in the limit Ug — ?■ oo, . Dotted line is Mi = 1 line, Eq. (A12); dashed line is the line 
of maximum deflection tan Xmax = l/2(tani?i) — cot i 
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In particular, for strong fluid shocks, o" = 0, in the hmit of high compression ratio, r/ ^ 1, 



^Uq sin ( 



1 sm 



tanx 



sin 2(j) 
2 + cos 2(1) 



(A13) 



Maximum deflection is reached at (/> = vr/S and is equal to x = "^/G (Konigl 1980). In the large 
compression limit the post-shock Mach number is 



Mf = 3 - (8/3)sin2</. 
The flow becomes sonic at the maximum turning angle (f) = tt/3. 



(A14) 



